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Abstract 

We explain how to apply Renormalization Group ideas to the analysis of the 
long-time asymptotics of solutions of partial differential equations. We illustrate 
the method on several examples of nonlinear parabolic equations. We discuss many 



applications, including the stability of profiles and fronts in the Ginzburg-Landau 



equation, anomalous scaling laws in reaction-diffusion equations, and the shape of 



a solution near a blow-up point. 
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1 Introduction. 



The development of a qualitative theory of infinite dimensional dynamical systems is 
a major scientific challenge. Such systems are expressed through (nonlinear) partial 
differential equations, and we shall concentrate on equations of the form 

m; = Am + F(m,V«,VVm). (1) 

where, u(x,t) G R (or C ), x £ R d , t G R + and F is nonlinear or random. Usually, the 
following problems are considered: 

1. Existence and regularity of the solution over a finite time interval. 

2. Extension of the first problem to an infinite time interval. 
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3. If the solution exists for all times, what is its asymptotic behaviour? How does the 
latter depend on the initial conditions? If the solution ceases to exist after a finite 
time, we may still ask about its behaviour near the point where it breaks down. 

The first question is often easy to deal with. In particular, perturbation methods 
(small initial data or small nonlinearity) tend to work well over finite (or short) time 
intervals. However, naive perturbation theory does not, in general, yield an answer 
to the last two questions. A similar situation was encountered in the theory of critical 
phenomena and in quantum field theory: perturbation theory works well whenever there 
are enough cutoff's, but tends to diverge when the latter are removed. In order to solve 
these problems, the Renormalization Group (RG) method was developed and proved to 
be very useful. One of our goals is to develop RG ideas in the theory of PDE's. Here, 
the role of the cutoff is played by the finite time. 

The RG method for nonlinear PDE's was developed by Goldenfeld, Oono and others 
i3|j . Previously, Barenblatt || emphasized the role of self-similar solutions of these 



equations. The latter can be viewed as fixed points of certain RG transformations (see 
below). Independently, we used RG methods to deal with an equation like (l),with 



a linear but random F [|| [L0[ |llf . Here we shall explain the basic RG ideas in the 
framework of PDE's, which turn out to be very natural and easy, and we shall review 
some of the rigorous results obtained so far. 

The general idea behind the RG approach is to solve the problem iteratively: first 
integrate the equations far away from the cutoff to be removed, i.e. here solve a finite 
time problem. Next, rescale the time and possibly the space and the u variables so as 
to produce a new problem similar to the original one. Finally, iterate and see what 
happens. The method will work when the new problem tends to be simpler than the 
original one (meaning, e.g., a weaker nonlinearity). Upon iteration, it will become 
simpler and simpler. Of course, this will happen only if we have chosen the scaling 
appropriately. 

The underlying properties of the dynamical system that make the RG method work 
are universality and scaling. In the PDE framework, scaling amounts to the observation 
that many problems have solutions that behave asymptotically as 

u(s,t) ~r*r$) (2) 

Often such a limit law is universal: the numbers a, (3 and the function /* will not depend 
on the initial conditions or even on the form of the equation. More precisely, pairs of 
initial data and equations will fall into universality classes, corresponding to given a, f3 
and /*. The fact that whole classes of problems may yield the same asymptotic behaviour 
is called universality. In fact, /* will be a fixed point of the RG transformation, while the 
exponents a and (3 are determined by the choice of the scaling. The basin of attraction 
of that fixed point is the universality class. The RG will lead to a dynamical system 
picture, with an analysis of the stable, neutral and unstable manifolds of the fixed points. 
The outline of the paper is as follows: first, we give an elementary definition of the 
RG transformation, starting from the (linear) heat equation. Then, we show how the 
RG flow acts on spaces of initial data and equations for various nonlinear perturbations 



of the heat equation, and for data that decay at infinity. In Section 3, we consider similar 
questions, but for data that do not decay at inifinity, which leads to the formation of 
patterns and fronts. Finally (Section 4), we apply the RG method to equations whose 
solutions blow-up in a finite time, and we give a list of possible profiles near a blow-up 
point. We end this review with a list of open problems. 



2 The Renormalization Group. 

2.1 The linear heat equation: Gaussian fixed points. 

To explain the RG idea, let us start with the simplest equation, the linear heat equation: 

Ut = u xx (1) 

in one dimension, with integrable initial data, u(x, 0) = f(x) G L 1 (R). The solution is 

u(x,t) = -L [ d yeM- { -^r-)f(y)dy (2) 

From this, it is easy to deduce that 

Vtu(tVi,t) - -=■ f(y)dy (3) 

V47T J 



as t — > oo (pointwise and in L 1 (R)). Let us reformulate this (trivial) result in the RG 
language. Notationally, it is useful to take the initial time equal to 1 instead of 0: 

u(x,l) = f(x) (4) 

This will not change anything to the long-time asymptotics. Next, take a number L > 1 
and define the renormalization map, acting on the initial data / as 

(R L f)(x)=Lu(Lx,L 2 ) (5) 

where u(-,L 2 ) is the solution of (1) at time L 2 . Thus, we solve the equation up to a 
finite time L 2 , rescale the solution and take it as a new initial data. Here, Rl is explicit: 
in Fourier transform, e tA multiplies by e~ tk and scaling x by L amounts to scaling k by 



-£. Thus, 



R L f(k) = e-^-^fiL-'k) (6) 



(we solve (1) over a time interval of length L 2 — 1] 
Observe that the set of multiples of 



f$(k) = e-« (7) 



forms a line of fixed points of the RG map (6). Note also that Rl satisfies the semigroup 
property: Since L is arbitrary, we may replace it by L n , and we have: 

R L n = R n L . (8) 

This is simply because, if we define Ul(x, t) = Lu(Lx, L 2 t), then 

(R L f)(x) = u L (x,l) (9) 

and ul satisfies (1) again. Equation (1) is scale invariant. 

Now, assume that our initial / belongs to the basin of attraction of the line of fixed 
points, i.e. 

Rlf - AfS (io) 

as n — > oo, for some A, in a suitable space. The basic observation is that this is 
equivalent to (3), by letting t = L 2n , A = /(0), and combining (5,8,10). Note that A 
depends on the initial data and that this dependence is the only trace left by u(x, 1) on 
the long-time asymptotics of the solution. This is why the latter is called universal. 

To see what kind of functions belong to the basin of attraction of this line of fixed 
points, write formally f(k) = /(0) + ck + • • •. Then, for k of order one, we have from 
the scaling 

Ihf(k) = f(0)f*(k) + O(L- 1 ) (11) 

and, for large k, we get contraction from the multiplication by e~ ~^' in (6). 
To be more precise, we use below the following norm (inspired by [p3[ ): 



= sup(l + |*|«)(|/(*)| + |/ / (*)|) (12) 

k 

on a space of functions with / e C 1 (R). We denote by B the corresponding Banach 
space. The value of q will be specified when we deal with nonlinear problems. For the 
moment, any q > will do. Any function / E B can be written in Fourier transform as 

f(k) = f(0)e- k2 +g(k) (13) 

with g(0) = 0. Then the convergence (10) with A = /(0) will follow from 

Lemma 1 There exists a constant C independent of L so that, if g G B and g(0) = 0, 

\\R L g\\<CL- l \\g\\ (14) 

Proof 

Write g{j) = Jo g'(p)dp so, \g(r)\ < VlMI- The scaling (6) brings a factor of L _1 in 
the derivative and \k\ n e~ k is bounded by a constant for any n > 0. □ 

Remark 1. At this point, we may comment on the choice of the constant L. It is 
arbitrary except that we take it large enough so that CL~ X < 1 in (14). This will be 
assumed, wherever needed, about any other numerical constant entering the proofs. We 
shall denote generic constants by C, even when they change in the same equation. The 
only important property of such constants is that they are independent of L. 
Remark 2. Using (13,14) we can prove (10), i.e. (3) for / 6 B. 



2.2 Nonlinear heat equations. 

Let us now consider a less trivial example, namely 

fh = u xx + \u p (15) 

where A could be scaled away but is introduced for convenience. Throughout this Sec- 
tion, we work in one dimension (to simplify the notation), but the extension to higher 
dimensions is straightforward. There is no problem in proving the existence of solutions 
for any fixed time, provided A is small enough, and we may define the renormalization 
group map as in (5). However, we do not have anymore the semigroup property (8), 
because the equation is not scale invariant. In fact, Ui{x,t) = Lu(Lx,L 2 t) satisfies 

u u = u Lxx + AL 3 "X (16) 

So, if we define Rl on the pair (/, A) by 

R L (f,X) = (u L (x,l),XL 3 - p ) (17) 

we recover the semigroup property (8), provided there is a set of pairs (/, A) that is 
mapped into itself by Rl. 

Now, consider p > 3. Then the line of fixed points (Af^O) of Rl is stable. So we 
may expect the long time behaviour of the solutions of (15) to be again given by Af£ 
for a suitable A. This is indeed the case. However, the constant A will not depend only 
on /, but also on p and A. 

We have the following result which is a variant of Theorem 1 in [[EJ. Fix p6N,j)>3 
in (15) and fix q > 1 in (12); let B\ be the unit ball in the corresponding Banach space 
B. 

Theorem 1 There exist an e > such that, if\\\ < e and f G B\, equation (15) with 
u(x, 1) = f(x) has a unique solution which satisfies, for some number A = A(f, X,p), 

)im\\Viu(Vi;t)-Af*(-)\\=0. 



Remark 1. For related results on the p > 3 case, see [|3], |35| (4^, |48| . 



Remark 2. The convergence in B implies convergence both in L and in L 



Let us explain the main ideas of the proof (for details, see ||12"|). The problem is to 
"construct" the constant A which is not given by an explicit formula as in (3). This 
construction is done inductively. 

First of all, the existence of the solution for finite times is rather trivial (we take A 
small and II f II < 1). Indeed, we can write 



u (t) = e^-^f + A /"* X dse sd2 u p (t - s 
Jo 



(leaving out the x dependence) and use a standard fixed point argument (see [|12j; here 
we need only q > 1, because there are no derivatives in the nonlinear term). Since the 
first term in (18) is the linear evolution, we may write 

Rf = R f + v (19) 

where R denotes the linear RG map (5) (we suppress here the L dependence of R) and 
||v|| is 0(A). Now, write 

f = A f* + g (20) 

with <7o(0) — 0. This splitting of / is not preserved by R, as it was in the linear case. 
Thus, we shall have to change A at each step. Write 

Rf = A 1 fS + g 1 (21) 

where again gi(0) = 0. From (19) and (20), we see that 

A = A + {)(0), (22) 

since R g (0) = (see (6)), and 

gi = Ro9o + v-v(0)f*. (23) 

Using Lemma 1 and \\v\\ = O(X), we have 

A 1 = A + 0(A) (24) 

and 

INI < CL- 1 ||( 7o ||+0(A) 

< L^ 1 ^ (25) 

for any 5 > 0, provided A is small enough, since \\g \\ = 0(1) for / e B\. 
Now iterate: under scaling A becomes Ai = AL 3 ~ P , and after n steps 

A„ = AL (3 " p)n (26) 

Thus, ||v n || will be 0(A n ) for the n th iteration, and, writing R n f = A n fQ + g n , we have 

A n+1 = A n + v n (0) (27) 

g n+1 = R° L g n + v n -v n (0)fo, (28) 



with the bounds: 



\A n+1 - A n \ < 0(A n ) (29) 

b„+i|| < CL- 1 ||( 7 „||+0(A n ) 

< L-^- 5)n (30) 
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(use (26) and p — 3 > 1 since p is an integer). 

Therefore there exists a constant A such that A„ — > A as n — > oo and 



\A-A n \<0{\ n ). (31) 



Then, 



i?7 = A n f*+g n 

= Af* + (A n -A)f* + g n (32) 

and 

\\Rlf - Afl\\ < L" Mn + 0{\ n ) (33) 

which goes to zero as n — ► oo. Using (5,8,17), this proves the claim, at least for a special 
sequence of times; the extension to all times is easy [JT2] . 

Remark 1. From (33,26), one gets an estimate 0(t~(~*~>) on the rate of convergence. 
Remark 2. The sign of A is not important here. It matters only for p < 3, see below. 
Remark 3. Rq is the linearization of the component of the operator R that acts on 
initial data, around the fixed point /g. All functions of the form k n e~ k ,n 6 N, are 
(Hermite) eigenvectors: 

R k n e~ k2 = L~ n k n e' k2 (34) 

The fact that the scaling transforms the linear evolution semigroup, e ta , which has 
continuous spectrum, into an operator with discrete spectrum, is one of the interesting 
features of this method. Using it, one can extend the above analysis and obtain higher 
order asymptotics of the solution in inverse powers of t, see subsection 2.5 below, and 



56, m 



2.3 The RG flow in spaces of equations and data. 

Of course, Theorem 1 can be extended to a much more general class of nonlinearities. 
Consider an equation of the form 

u t = u xx + F(u, u x , u xx ) (35) 

Then, the map Rl will be defined in general by 

R L = (L a u L (-A),F L ) (36) 

where 

F L (a, b, c) = L 2+a F{L~ a a, L^^b, L~ 2 ' a c) (37) 

If we consider a monomial F(a,b,c) = a Ul b n2 c n3 , and the previous scaling a = 1, we 
have F L = L~ dp F with 

dp — n\ + 2n 2 + 3n 3 — 3 (38) 



and using standard RG terminology, we call the monomial relevant if dp < 0, marginal 
if dp = 0, and irrelevant if dp > 0. 

Considering n, > 1 or n, = 0, the nonlinearities are easily classified: u p is relevant 
for p < 3; w 3 and Mtt x are marginal, and everything else is irrelevant. Then, Theorem 1 
extends to the irrelevant case, provided we take q > 3 in (12) and F : C 3 — > C analytic 
in a neighbourhood of zero [0 . 

Remark. More general irrelevant nonlinearities (nonlocal, integro-differential terms) 
were considered by Taskinen [53]. Besides, as mentioned in [T2|], we may replace the 



second derivative in (15) by —(—A)' 3 / 2 and obtain similar results. 

A nice extension to the study of waves satisfying a generalized KdV-Burgers equation 
was made in |J. There, the highest order (third) derivative turns out to be irrelevant. 

The two marginal terms exhibit different behaviours. For w 3 , it is essential to replace 
the plus sign by a minus one (with A > 0) in equation (15) (with a plus sign, and p < 3, 
any positive initial data leads to a solution that belows up in a finite time |]3"4"|, |55|). The 



nonlinearity turns out to be irrelevant in the next order. The coefficient A n satisfies a 
recursion of the form A n+ \ = A n — A/3 A 3 + higher order terms, where (3 > is an explicit 
constant. This implies that A n ~ n~z (here we use A positive) and since t = L 2n , this 
translates into the following logarithmic correction (see |I2"|): 



u 



(£VM) ^ (tlog*)-'/o(0 (39) 



On the other hand, for uu x , we have Burgers' equation and, via the Cole-Hopf trans- 
formation, we get a new line of fixed points: 

r( x ) = A-\og(l + Ae(x)) (40) 

where e(x) = J**^ e 4 dy is the error function and A is a parameter. A theorem analo- 



gous to Theorem 1 can be proven for these new fixed points [12|. Of course we can also 



add to these marginal terms a general irrelevant F and obtain a similar result. 

2.4 Non- Gaussian fixed points. 

Let us now turn to the relevant case, u p , 1 < p < 3. Actually, since we will not restrict 
ourselves to integer p or to everywhere positive w's, we shall consider the equation: 

u t = u xx - w|w| p "\ (41) 

but the minus sign will again be essential. Here, the parameter A of (15) has been 
eliminated by a rescaling of u. 

Since the nonlinear term is relevant, the Gaussian fixed point will be of no use. In 
order to find another fixed point, consider the transformation (36) with a = -^-r. The 
reason for this choice is that equation (41) is then scale invariant: Fp = F. It is easy 



to see that we shall have a fixed point of the transformation if we have a solution of the 
form: 

u(x,t)=t~p=if(xt~*) (42) 

where /(£) is a function of one variable which solves: 

/" + Uf + ^—-f p = 0. (43) 

2 p — 1 

These solutions are called self-similar. The theory of positive solutions of (43) has been 
developed in R |35], JS|. The main result is that, for any p > 1, there exist smooth, 
everywhere positive solutions, / 7 , of (43) with / 7 (0) = and / 7 (0) = 7 for 7 larger than 
a certain critical value 7 P (but not too large). The decay at infinity of these solutions is 
given by 

f\(x) ~ |x|~5=i (44) 

as \x\ — > 00 if 7 > 7 p , while, for 7 = 7 P , it decays at infinity as 

f lp (x) ~ Ixl^ 71 " - e~~. (45) 



Remark. The existence of a critical 7 P can be understood intuitively by viewing (43) as 
Newton's equation for a particle of mass one, whose "position" as a function of "time" 
is /(£). The potential is then U(f) = 2 < -v> ~~ ^~+T anc ^ ^ ne "friction term" |/ depends 
on the "time" ^. Hence, if / 7 (0) = and / 7 (0) = 7 is large enough, the time it takes to 
approach zero is long and, by then, the friction term has become sufficiently strong to 
prevent "overshooting" . However, as p increases, the potential becomes flatter and one 
therefore expects 7 P to decrease with p. 

We shall now explain why these self-similar solutions are stable. Consider the initial 
data (taken as before at time 1 for convenience) 

u(x, 1) = / 7 (x) + h(x) (46) 

with 7 > 7 P and h G B, where B is the Banach space of L°° functions equiped with the 
norm (with some abuse of notation!) 

11/zlU = esssup|/z(^)(l + |e| g )|. (47) 



One can show the following p6 : 



Theorem 2 Let 1 < p < 3. There exist e > 0, C < 00 and \x > such that, if the initial 
data u(x, 1) of (41) is given by (46) with h e B and satisfies 

, < e, or h > 0(a.e.) 



then, (41) has a unique classical solution and, for all t, 

||t^ u (.t3,t)-/ 7 (.)||oo<cr'*||/i| 



Remark. For related results on the stability of self-similar solutions, see E9L BQ, 31, 35 
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Let us explain the main ideas of the proof. Given the initial data (46), it is convenient 
to rewrite (41) in terms of the variables £ = xt~^ and r = logt; so, define i>(£, r) by: 

u(x,t) = t~^ =l {f 1 {xt~^) + v(xt~2 ,logt)) (48) 

where now 

«(e,0) = /i(O- (49) 

Then, (41) is equivalent to the equation 

v T = Cv - (|/ 7 + w|"- 1 (/ y + v) - / 7 P - p/*" 1 *) ee Cv + N(v) (50) 

where we used the fact that / 7 solves (43) and we gathered the linear terms in 

with 

d 2 £ d 1 

c ° = ie + U + —i' (51) 

and 

To prove the theorem, it suffices to show that the corresponding solution of (50) goes 
to zero as r — > oo. For that, the main estimate will be that the semigroup e rC contracts, 
at least for r large: 

||e r£ || < Ce-» T (53) 

for some \x > 0, C < oo. Given (53), the control over the nonlinear term in (50) is 
standard. 

There are two important ingredients in the proof of (53). The first is the fact that 
e rC is a contraction in a suitable Hilbert space of rapidly decreasing functions. To see 
this, note first that £ can be conjugated to the following Schrodinger operator: 

?r -£ d2 ¥ 1 l 

e8£oe '=dp-i6-i + ^r (54) 

10 



i.e. to minus the Hamiltonian of the harmonic oscillator. Thus, Co is self-adjoint on its 
domain T>(C ) C L 2 (R,dfi), where 

i. 
d/i(C) = e*d£. 

Co has a pure point spectrum {— j — ~ — — | m = 0, 1, • • •} and the largest eigenvalue 
-^- — | is positive if 1 < p < 3. Thus e rCo is not contractive and, for e r( - Co+v ^ to contract, 
we need to use the potential in a non-trivial way (this is the reason why 1 < p < 3 is 
harder than the p > 3 case). 

Remarkably, it is possible to prove that C < —E < without a detailed study of the 

function / 7 , but only using equation (43). Indeed, one first shows that C is self-adjoint 

and that — E 7 , its largest eigenvalue, satisfies — E 1 < —E lv . So, writing E = E 7p , it is 

enough to show that —E < 0. Next, it is easy to see, using the Feynman-Kac formula 

61] and the the Perron- Frobenius theorem j|T) that C has a unique eigenvector Q with 



eigenvalue —E and tt can be chosen to be strictly positive. Write (43) as 

Cf = -(p-l)f p . 
So, 

(n,£/) = -(p-i)(n,/*) 

where (•, •) denotes the scalar product in L 2 (R, dp,). By the self-adjointness of C and 
the definition of Q, i.e. CQ = —EQ, we have 

since f2 and / are strictly positive. This proves our claim. 

Notice that functions in L 2 (R,, dp) have essentially a Gaussian decay at infinity, which 
is much faster than what is allowed in our Banach space B, see (47). An extra work is 
needed, using the Mehler's formula for the kernel of e TjZ ° (see |16J for details). 

Finally, let us express this result in the RG language. The RG map has a fixed point 
i2i(/ 7 , F*) = (/ 7 , F*) for F*(u) = — it|it| p_1 . There are different fixed points for different 
values of 7, but since the decay of the functions in B is faster than the one of any of 
these fixed points, each of them is unique in the corresponding set {/ 7 + h} for h as in 
Theorem 2. In other words, unlike the situation for p > 3, we do not have to deal with 
a line of fixed points and here there is no constant like A in Theorem 1. This is reflected 
in the fact, (53), that the linear semigroup e rC contracts. 

Remark. The reader who is familiar with the theory of critical phenomena will recognize 
an analogy between the behaviour of an Ising model at its critical point, when the 
dimension of the lattice changes and the asymptotic behaviour of the solution of (41) 
when p varies. This is further discussed in ||12|| . Here, our results for p < 3 are non- 
perturbative, i.e. we do not use any "e-expansion" (see |6f| for an analysis of the 
bifurcation at p = 3). 

Also, the reader may notice that our approach to the renormalization group is close to 
the "Wilsonian" one in the theory of critical phenomena, while the method developed in 

11 



43| is analogous to the renormalized perturbation of quantum field theory. For a rigorous 



treatment of Barenblatt's equation, which was one of the first example analyzed in [43 

sec 



2.5 Higher-order asymptotics. 

It is interesting to reexamine the p > 3 case, using what we learned from the analysis of 
p < 3. This has been done by Wayne |66 |. Namely, consider equation (15), and make 



the change of variable (48), but replacing / 7 by 0, i.e. 

u(x, t) = t~~v(xt~ 1/2 , logt) (55) 

Then t>(£,r) satisfies 

v T = C v + \v p (56) 

with Cq as in (51); we have computed the spectrum of £ which is {~zj — \ 



2 



\m 



0, 1, • • •}. Now, notice that, for p > 3, this spectrum is entirely negative. Thus e TC ° 
contracts exponentially and it is not hard to show that, for A small, and v (£, 0) bounded 
in a suitable norm, the same contraction holds for the solution of (56). Actually, since 

the largest eigenvalue of Cq corresponds to m = 0, we expect a decay like e ^p- 1 2> , 
which, using r = logt and combining it with the t - " prefactor in (55), leads to a £ -1 / 2 
decay of u(x, t), for all p > 3. This of course has to be the case, in view of Theorem 1. 

Actually, much more can be said. Take the n largest eigenvalues {Ai,---,A n } of 
Co and the corresponding subspaces. Then, using a theorem of Gallay |J7), Wayne 
constructs an invariant manifold (in the v variable) such that any solution (for p > 3 
and small initial data) will converge to that invariant manifold at a rate at least of order 
-/--(Xn+i-S) (-^i-th, as before, S arbitrarily small for sufficiently small initial conditions). 

Since, for n = 1, the eigenvector of Cq is e~^ / 4 , one easily recovers Theorem 1. But 
all higher order asymptotics in time of v (£, r), hence of u(x, t), can also be obtained in 
terms of Hermite functions (such higher order corrections were also derived in |5q| ). 

For p < 3, an invariant manifold can still be constructed but it becomes unstable 
(— |-j- — | > 0) and thus, it does not, in general, give information on the long-time 
asymptotics. Of course, we know from Theorem 2 what the situation is: instead of 0, 
we have to consider the solution / 7 and, presumably, a similar picture (with invariant 
manifolds and higher order asymptotics) holds there, since the spectrum of the operators 

C for p < 3 is qualitatively similar to the one of Cq for p > 3. 

2.6 An application to reaction-diffusion equations. 

In [f|] , Berlyand and Xin consider the following model which leads to a nice RG analysis 
and which exhibits the novel feature of anomalous exponents depending on the initial 

data: 

u t = Au + vuP- 1 (57) 

v t = A.~ 1 Av-vu p ~ 1 (58) 

12 



where u = u(x,t),v = v(x,t),x G R n (we take n — 1 for simplicity) and A > is the 
Lewis number. Equations (57,58) model a chemical reaction A — »■ B where v is the mass 
fraction of reactant A and u the one of reactant B. 

First, observe that for p > 3, the analysis done for (15) applies and the asymptotics 
of u,v is Gaussian. So let p = 3, so that the nonlinear terms are marginal, and let us 
see heuristically what happens. If u and v are positive (which is physically necessary 
and is preserved by (57,58) due to the maximum principle) then, again by the maximum 
principle, u is larger than the solution of the heat equation with the same initial data: 

u(x,t)>^e-£. (59) 

where A depends on u(x, 0). 

Then, v is less than the solution of (58) with w p_1 = u 2 replaced by its lower bound 
(59). Calling v this upper bound on v, v solves a linear equation 

v t = A^Av - AH^e'^v. (60) 

It is not hard to see, by direct substitution, that (60) admits a self-similar solution of 
the form 

v{ X) t) = t^a^). (6i) 



Vi J 



Indeed, writing /*(£) = e s *!>(£), \I> solves 



- A" 1 *" + de + fe-^ 2 + 1)* = f * (62) 

16 An 4 2 

We see that, requiring /* to be positive means that \l/ is the ground state of the operator 
in the LHS of (62) i.e. of an harmonic oscillator perturbed by a positive, rapidly decaying, 
potential. The exponent ^ is the corresponding ground state energy. For the harmonic 
oscillator, \ — \ and, perturbation theory tells us that, for A small, 

a 1 A 2 , A 1 

- = - + + h.o.t. in A > -. (63) 

2 2 47TV2A- 1 + 1 2 v ; 

Thus, a depends on A, i.e. on u(x, 0). Since a > 1, v, and hence v, decay strictly faster 
than the solution of the heat equation. Inserting this in (57) gives an upper solution for 
u, which solves an equation where the nonlinear term is now irrelevant. 

These considerations lead us to expect u to have the heat equation decay: 

u(x,t) ~ A e - 2 /4< ( 64 ) 

and v to have the anomalous decay: 
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where A, B and a depend on the initial data. Upper and lower bounds of the form 
(64,65) are proven in ||. From the RG point of view, one can define a map 

where 

Ul(x, t) = Lu(Lx,L 2 t), (66) 

v L (x,t) = L a v(Lx,L 2 t). (67) 

We have thus a two-parameter family of fixed points (AfQ,Bf*, A j), and the scaling 
exponent varies continuously with A. 

3 Patterns and fronts. 



3.1 The Ginzburg-Landau equation. 

In the previous section, we have seen that nonlinear equations with initial data decaying 
at infinity produce universal, Gaussian or non-Gaussian, diffusive profiles. Here we shall 
show that other types of asymptotic behaviour can also exhibit such universality. To 
discuss a concrete example, consider the Ginzburg-Landau equation 

u t = u xx + u-\u\ 2 u (1) 

where u : Rx R — ► C is complex. This equation has a two parameter family of stationary 
solutions 



u qe (x) = sjl-q 2 e i ^ + ^ (2) 

and a natural question is to inquire about the time development of initial data u(x) 
which approach two such solutions at ±oo: 

Jirn^ \u(x) - u q±e± {x)\ = 0. (3) 

This problem has been extensively studied for equation (1) with u real and u{— oo) = 
1, u(oo) = (i.e. q- = 0, q + = 0) |T], 0. Then the solution takes the form of a 
propagating front. This occurs because u(x) = 1 is a stable stationary solution, while 
u(x) = is an unstable one. For complex u, the solutions (2) are stable, under small 
perturbations, for q 2 < | (the Eckhaus stable domain) |23|| . 

We shall consider two questions. The first one was suggested in ||24|| , namely we 
take q± in (3) small, belonging to the Eckhaus stable domain, but not necessarily equal. 
What is then the long-time asymptotics of the solution? The second question concerns 
the stability of the real front solutions of (1) for real or complex u. For reviews on these 



questions, we refer the reader to p2|, g5 . 
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3.2 Non- Gaussian patterns. 



For the first question, we considered in |13] a class of initial data satisfying (3) and we 
showed that, for any interval J, 



C 

sup \u(x,t) - e lv/ *^*M g »6)*(x)| < — j= 
xei yt 



(4) 



where the constants q*, <fi* and 9* depend only on the boundary conditions (3). For a 
more detailed bound, see Section 3 of |Tj| . Graphics of the solution can be found in |25| . 



In order to see where (4) comes from, we write, following |24 

u= (l-s)e^. 
Equation (1) becomes, in these variables, 

St = s xx -2s + 3s 2 - s 3 + 4> 2 X - s<f) 2 x = s xx -2s + F(s, Y . f . , 



\s,(p x ) 



^(pxS x 



= (f>xx + G(s, s a 



(5) 



(6) 



with the initial data ( taking again t — 1 as initial time) 

lim s(x, 1) = s± , lim \<p(x, 1) — (f)±x — 9±\ =0. 

x— >±oo x— >±oo 

where 2s± = F(s±, <p±). We proved in |^3| the following asymptotics, as t — > oo: 



(7) 



s(x,t) 



^ + ^f'^ + < ] 



(9) 



for a set of initial data satisfying (7), with (f>±, and 8± small enough (see again [13 
for the precise definition of the norms in which these limits take place). The functions 
0*, 6*, s*, r* are universal, depending on the initial data only through the boundary 
conditions (7). They are smooth and therefore (expanding 0* to first order, and using 
(2)) the phase of u will have the asymptotics (4), with 

(f)* = <f(0), q* = <j)*'(0), 6* = 0*(O). 

The behaviour of s will be discussed below. 

The peculiar scaling behaviour exhibited by (8) can already be understood through 
the linear heat equation. Previously, we considered the Gaussian fixed point of the RG 
in a space of integrable initial data. Since the heat equation has stationary solutions 
of the form a + bx, we may also consider the solution of <\>t = <p xx with initial data 
4>(x, 1) = f(x) satisfying (7). The solution can be computed: 



r^(j)(Vix,t+ 1) 
-X + (6 



'A<Kt 



,-\{y-*) 2 



> T '" "' f(Viy)dy, 



t— »oo 



-^){x + 2—)e{x)=4>l{x) 



(10) 
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where e(x) = /f^ e~~dy is the error function. This gives the first term in (8), with 0* 
replaced by 0q. If we want to obtain the next correction, we write 0* = 0q + 9, where 
9(x,t) solves 9 t = 9 XX with 9(x, 1) — ► 9± as x — ► ±oo. One computes that 

fl(x,t) - 91 [^ = 9. + (9 + - 9_)e (^ (11) 

and this yields the second term in (8). It is trivial to check that 0q (resp. 9q) is a fixed 
point of the RG map (2.36) with a = — 1 (resp. a = 0), and F — 0. 

We shall explain below that 0*, 9* are small perturbations of 0q, #o- However, as we 
saw, the scaling (8,4) can be understood qualitatively on the basis of the heat equation, 
and the boundary conditions (7). We still need to understand (9). 

The latter comes from the "slaving" of s to 0, due to the linear —2s term in (6). 
This will imply that s will have the right form so that (4) holds. Concretely, the slaving 
means that, to leading order, we may solve the algebraic equation —2s + F(s, <f) x ) = 0, 



which is equivalent to 1 — s = Jl — (4> x ) 2 , and substitute the result in the equation (6) 
for 0. Since the derivative of s is proportional to the second derivative of 0, we obtain 
an equation for only, of the form: 

t = (1 + a{4> x ,(t> xx ))(t> xx , 0(x, l) = f(x) (12) 

with the boundary conditions (7) for the initial data /. The function a is analytic around 
the origin. 

The RG map acts again on pairs of initial data and equations and can be defined 

(with a = — 1) as 



where 



and 



R L (<f>,a) = ((f)L,a L ) (13) 



<f> L (x,t) = L~ L (f)(Lx : L 2 t) 



cll(u,v) = a(u, L 1 f). (14) 

The semigroup property follows by observing that 0l satisfies 

4>Lt = (1 + a L ((f) Lx , 4>Lxx))<t>Lxx- (15) 

We may now iterate R as before, i.e. solve a finite time problem, to study the 
asymptotics of (12). One wants to show that there exist functions 0*,a* such that 

it>2(0,a)^(0*,a*) (16) 

where 

i? L (0*,a*) = (0*,a*) (17) 
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is the fixed point of the RG, corresponding to the scale-invariant equation <f> t = (1 + 
a>*) ( l ) xx- Then, replacing A= by £, the asymptotics of the original problem is given by 

0(^ 1/2 ,*)~^0*(O- (18) 

Because of the factor L _1 in the second argument of a in (14), the fixed point equation 



LS 



<j>t = {l + a*{<f> x ))<i> xx (19) 

where o*(-) = a(-,0). The fixed point is the scale invariant solution <fr(x,t) = \/t(f)*(A=). 
We get for 0* the equation 

with a* = a* (4:0*) and, for small <fi±, we look for a solution 

0* = 0S + P (21) 

where p(±oo) = and 0q is the "Gaussian" solution (10), which solves (20) with a* = 0. 
This is easy to solve by a fixed point argument (see Proposition 1 in [[TI| or the Proposition 
in Section 4 of [|T2l ). 



This gives us the first term in (8). Turning to 8*, we write 

<P(x,t) =<f)*{x,t) + 6{x,t) (22) 

where, with an abuse of notation, <p*(x,t) = \^t(p*(A=) and 0* is given by (21), while 
solves (12). Then, 

Qt = Qxx + (a(j) xx - a*(fr* xx ) (23) 

with #(±oo) = 9±. Now we set 

6 L (x',t')=6(Lx' J L 2 t'). (24) 

Putting x = Lx',t = LH\ we have J^0*(Lx', L 2 t') = L l - l -^(jf(x' , t') and therefore 9 L , 
satisfies the equation (replacing (x',t') by (x,t)): 

0u = B Lxx ^L{a^l^L-H Lx ,L-^l x ^L- 2 Q Lxx W xx ^L- x B Lxx ) 

-o*(C)«- (25) 

Thus, reasoning as above, we expect 

6 L n -► 6* (26) 

where 9*(x,t) = d*(-j=) satisfies the L — > oo form of (25): 

0t = 8*xx + a *°*xx + M<P* X , 0)0* + a v ((f)* x , 0)0^)0^. (27) 
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This is a linear equation, easy to solve, whose solution is, for 9± small, a small pertur- 
bation of the "Gaussian" solution (11) (which solves (27) with a = 0). Finally, one sets 
(with the same abuse of notation) 

<f)(x, t) = <j>*(x, t) + 0*(x, t) + V(x, t) = Vi<j>*(-^) + 0*(-^=) + *l>(x, t) (28) 

As for the s variable, one gets 

s(x, t) = s*(x, t) + r*(x, t) + v(x, t) (29) 

where 

s*( x ,t) = s*(^=), r *(x,t) = ±=r*(^=) (30) 

are fixed points "slaved" respectively to <p* and 9*. They satisfy the boundary conditions 

lim |s*(x) — s+\ = , lim r*(x) = . (31) 

The equations for ip and i> are rather complicated, but are essentially of the form 
heat equation plus irrelevant terms, in the sense of the previous section. Since now ip 
and v decay to zero at infinity, we are in the situation of that section and, taking small 
initial data (in a suitable norm), one shows that the corresponding solution diffuses to 
zero. This, then, allows us to prove equation (4). 

3.3 Stability of Fronts in the Ginzburg-Landau Equation. 

Let us write the Ginzburg-Landau equation (1) in radial and angle variables, u = re l ^\ 

r t = r xx + r(l-<f>l)-r 3 (32) 

(fit = <\>xx + 2r~ 1 (j> x r x (33) 

and let us discuss the stability of the front solutions of these equations. It is well known 
that these equations have real, positive, front solutions, i.e. solutions of the form 

= O,r = r c {x-ct) > 0, (34) 

such that r c interpolates between a stable and an unstable solution of (32), i.e. r c — ■> +1 
for x — » — oo,r c — > 0, for x — ► +oo. Indeed, from (32,34), we see that r c satisfies 

r c + cr c + r c — r c = (35) 

which, if we reinterpret the variable as "time", can be seen as Newton's equation of 
motion of a particle of mass one subjected to a friction term cr' and to a force deriving 

2 4 

from the potential y ~~ x> w hich is an inverted double- well. It is intuitively clear and 
easily proved that, for c not too small, solutions exist that satisfy the required conditions, 



i.e. such that r c tends, as "time" goes to +00, to zero, the stable critical point of the 
potential, and to one as "time" goes to —00. For large "time" u = x — ct,r c (u) will 
decay exponentially, as is seen from the linearization of (35) at r = 0. One gets 

r c (u) < (d + C 2 u)e^ u (36) 

where 7 is given by 7 2 — cy + 1 = i.e. 

1 



7c = -(c-v / c^4) (37) 

which is real for c > 2, in which case 7 C < 1 (actually, one can take C 2 = in (36), 
if 7 < 1). Thus, the larger the friction, the slower the decay. For c < 2, the solution 
"overshoots" the minimum at zero, i.e. r c is no longer positive. Each of the solutions 
r c with c > 2 is stable under real perturbations (0 = 0): if we start with initial data 
r(x, 0), with r = r c + s with < r < 1, s decaying faster than e~ lcX for x — > +00, r(x, t) 
will converge, as t — > +00, to r c (x — ct), see |1|, ||, |22| . 

However the solution with c = 2, 7 C = 1 is more stable than the others in the sense 
that any initial data r(x,0) with < r < 1 which decays faster than e~ x as x — > +00 
(in particular, if r is of compact support) will converge, as t — ► +00, to r 2 (x — 2t) |l|, ||. 
From now on, we shall concentrate on the most interesting front, namely the one with 
c = 2, 7 = 1 and we write r for r 2 . 

We consider a complex perturbation of r : r(x, 0) = r + s with 0(x, 0) 7^ and 
<p(x, 0), s(x, 0) are small in a suitable sense. The equations satisfied by (ft and s are : 

4>t = <\>xx + 2{r x <f> x + s x (j) x ) (r + s) _1 (38) 



s t = s xx + s(l - 3r 2 ) - rcg - s(g - 3rs 2 - s 3 (39) 

Since r is a function of x— It it is convenient to consider also the equation in the frame 
of reference of the front : let u = x — 2t and 0/(w, t) = (p{u + 2t, t),Sf(u, t) = s(u + 2t, £); 
then (fif and Sf satisfy equations like (38,39), with 2(pf u added to the RHS of (38) and 
2sf u added to the one of (39). Now r = r{u) is time-independent. 

To understand the expected behaviour of 4>f(u, t), let us consider the linearized equa- 
tion around the zero solution : 

<t>ft = <Pfuu + 2r u (f> fu r~ 1 + 2(j) fu . (40) 

It is convenient to rewrite this equation as a heat equation with a potential: Let 

<j) f (u, t) = e-"r(u)-VK *)• ( 41 ) 

Then, ip satisfies 

A = ^uu - Vi> (42) 
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with 






V = 1 + — + 2- = r\ (43) 

To derive the last equality, we used eq.(35), satisfied by r. Since r ~ 1 for u — > — oo, 
r ~ e~" for u —>■ +oo, we have V ~ 1 for u — > — oo, V ~ for w — > +oo. 
So, starting with ip(u, 0) localised around u = 0, we expect that 



* « 2 



U — >• — OO 



V>M)~<( $ (44) 

w — > +oo. 



Vt 
Hence, using (41) and the asymptotic behaviour of r(u 



0/(M)~ f " ( 45 ) 



-t—^r-u 
e « 

u 

e ~%t 

Vt 


u - 
u - 


-> — oo 
-> +oo 


(u+2t) 2 
e it 

ft 


u - 


-> — oo 


e ~%t 


u - 


■*■ +oo. 



which can be written as 

<mm)~<! y* '^ ^ (46) 

Since u + 2t = x, the first part of (46) is a diffusive wave stationary in the fixed frame. 
The second part represents a diffusive wave which is "carried along" by the front. This 
is a rough, but basically correct picture. 

Let us consider the linear equation for Sf, in the front frame : 

Sft = s fuu + s fu + s f (l-3r 2 ) (47) 

Writing Sf = e~"cr, we get 

°t = Ouu - Va (48) 

with V = 3r 2 . Following the analysis leading to (44), we get 

,,2 
s -3t-^-u _ £ -2t 



e « as u — »■ — oo 



Vt V~t 

e~ u as u — > +oo 



s(u,t)~{ J 1 ^ (49) 

e "IT 

Vt 



There is a "wave" which is stationary in the front frame, but exponentially decreasing 
in u , while the wave which stays in the fixed frame is suppressed by the factor e~ 2t . 

The rigorous results that justify this picture are of two types: Using the RG method, 
Gallay |3J| was able to obtain very precise asymptotics on how a small perturbation of 



real fronts diffuses to zero (this improves previous results of |53], |6T|). So Gallay considers 
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equation (39) with <fi = 0, in the front frame, whose linearization is given by (47). He 
writes 

s(u,t) = r'(u)w(u,t) (50) 

and studies the behaviour of w(u,t) as t — ► oo. Since r(u) goes exponentially to 1 or 
as u goes to — oo or +oo, r'(u) will be localized around u = 0. The main result is that 
w has the following universal asymptotics [Q: let u = £,\/i, then 

w(tVi,t)~Ar 3 / 2 r(0 (5i) 

where 

/ (c)-| e - € V4 ife>0 

and A depends again on the initial conditions. The limit (51) holds in a weighted L l r\L°° 
norm. Going back to (47-49), the power £~ 3 / 2 is easy to grasp intuitively: the potential 
V in (48) plays the role of a barrier around u = (V goes exponentially to 3 for u — > — oo 
and exponentially to far u — > +oo). The simplest approximation is to replace the RHS 
of (48) by a Laplacian on R + with Dirichlet boundary conditions at u — 0. And that 
operator leads to a £~ 3 / 2 decay of the solution. This effect was not taken into account 
in (49). 

The second type of results deals with the complex perturbations. In [H j , we consider 
initial data in the Banach space of C 1 -functions <ft, s with the norm 

OH = sup(l+ | x |) 3+5 (| <f>(x) I + I 4/{x) I +(1 + e*)(| s(x) | + | s'(x) |)) (52) 



and prove 

Theorem 3. For any 5 > there exists an e > such that equations (38, 39) with initial 
data <p(x,l) = <p(x),s(x,l) = s(x), and ||(0,s)|| < e, have a unique classical solution 
4>(x, t), s(x, t), for all t > 1, such that 

\4>(x,t)\ < H+* (53) 

(l + e u )\s(x,t)\ < r 1+5 (54) 



Remark 1. The power laws of the decay in time are presumably optimal (except for the 
5) and are different from those of Gallay, because the diffusive wave that is stationary 
in the fixed frame, for the variable (see (45)), goes only diffusively to zero. This in 
turn slows downs the decay of the s variable, due to the nonlinear term s<p 2 x in (39). 
Remark 2. The nonlinear terms in (38,39) turn out to be irrelevant, in the RG sense. 
However, this is not a simple affair: to show it, one has to take into account the precise 
decay of (f) and s both in the fixed and the front frames. This makes the proofs rather 
complicated. 

Remark 3. Finally, let us mention that Eckmann and Wayne [|||, using a completely 
different (and simpler) method, namely coercive functionals, have proven similar results: 
they can consider a larger space of perturbations (s, (f>) than the one defined by (52), 
but they do not obtain explicit upper bounds on the decay in time. 
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4 Universality in Blow-Up for Nonlinear Heat Equa- 
tions 

4.1 Statement of the problem. 

Let us now consider equations for which global existence results do not hold: the solu- 
tions of the initial value problem 

u t = u xx + u p (1) 

where p > l,u = u(x,t),x G R, and u(-,0) = w G C°(R), will, for a large class of 
initial data uq, diverge in a finite time at a single point (for reviews on this problem, 
see [0, ^5f). Again, we limit ourselves to one space dimension, but the generalization 
is straightforward. 

The RG ideas can be applied to the analysis of the profile of the solution at the time 
of blow-up. To explain what this means, let us fix the blow-up point to be and the 
blow-up time to be T. Then, we ask whether it is possible to find a function f*(x) and 
a rescaling g(t, T) so that 

lim(T -t)^ u(g(t, T)x, t) = f*(x) (2) 

Moreover, we want to see how g or /* depend on the initial data. 

The prefactor (T — £)" in (2) can be understood easily: for initial data uq(x) 

constant in x, u(t) solves the ODE u t = u p , i.e. u(t) = ((p — 1)(T — t))~ for T = 
(p — 1) -1 'U ~ P . In [|4], [45], [|6|, |64| (see also |32[ |33|) several possible /*'s are discussed, 



and the set of initial data that will lead to a given /* is partially characterized. In |15 
we showed that there exists, in the space of initial data C°(R), sets Aik of codimension 
2k, such that, for uq G A4k, the limiting behaviour (2) is obtained, in the case k — 1, for 

g(t,T) = ((T-*)|log(T-*)|)' (3) 

f*(x) = (p-l + Fx 2 ) 1 ^ (4) 

where b* = , - , and in the case k > 1 for 

g(t,T) = (T-t)^ (5) 

f* b (x) = fp-l+bx 2 ")^. (6) 



where now b is an arbitary positive number. As shown in [15], one can also add suitable 



(irrelevant) terms to (1) without affecting the result. It was shown in [[44], |46|, g4| |33[ 



that, under quite general hypotheses, (3,4) or (5,6) are the only possibilities (see also 



[65] for a formal analysis). Moreover, solutions that behave as in (5,6) for k = 2 were 



constructed in |H4 



The codimension of M.k for k = 1 is easy to understand: since we have fixed the 
blow-up point (to zero) and the blow-up time (to T), we have to fix two parameters in 
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the initial data. To reach the other profiles, 2k — 2 additional parameters need to be 
fixed in the initial data. The k = 1 situation is therefore the most generic one. 

In the RG language, /* and /£ can be viewed as fixed points of a renormalization 
group transformation having 2k unstable ("relevant", in renormalization group termi- 
nology) directions. Thus, to converge towards the fixed point, one has to fine-tune 2k 
parameters (one for each unstable direction) and this explains why M.^ is of codimension 
2k, and in what sense /*, f£ are "universal". In addition, we encounter also one neutral 
("marginal") mode, which, for k — 1, turns out to be stable when nonlinear effects are 
taken into account and for k > 1 parametrizes a curve of fixed points. 

Our results are perturbative, i.e. the sets M.k consist of initial data that are close to 
the corresponding fixed point. Therefore, our results are similar to those of Bressan || 
who considers a nonlinearity e u instead of u p and obtains the universal profile analogous 
to our k — 1 case. The connection of blow-up and center manifold theory was used 
earlier in the work of Filippas, Kohn and Liu [B2, [33[ and of Herrero, Velazquez [HI, HR 



46, 47, 64] and Galaktionov |6q|. Futhermore, the scaling and the dynamical systems 



aspect of our work goes back to Giga and Kohn |38|, |39], |40| . Rescalings as in (7) below 



were used as a technique for numerical computation in ||. 

Let us first describe a change of variables that transforms the problem (1) into a 
problem of long time asymptotics: we write (1) in the "blow-up-variables": given a 
m:Rx[0,T)^R, define : R x [- log T, oo) -> R by 

u(x, t) = (T- t)-^cj>( * - log (T - £)). (7) 

Then u is a classical solution of (1) if and only if 0(£, r) is a classical solution of 



-2 



1 . 1 



L-'fee - :tt£0£ -<i> + 



2A; ? p-l 
0(£,T O ) = T^o(T^) (9) 

where r = — logT, and 

L T = eH 1 " 1 /^ (10) 



We construct in ]TJ| global solutions of (8), with suitable initial data, thereby estab- 
lishing blow-up for (1). Note that, for k — 1, the scaling in (7) differs from the one used 
in (3) by a factor | log(T — t)] 1 ^ 2 = r 1 / 2 . Actually, the situation where k > 1 is easier to 
understand heuristically, so let us start by discussing the latter. 

4.2 Analysis of k > 1. 

For k > 1, as r — > oo, the factor L~ 2 in front of the second derivative in (8) leads us to 
consider the solutions of 

^"^-^T^ (u) 
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Observe that the "fixed points" f£(x) = (p — 1 + bx 2k ) 1 p of (6) are stationary solutions 
of that equation. The latter can of course be integrated in closed form, but before doing 
that, let us first look at its linearization around the constant solution <ft = (p — 1)~. 
The linear problem is (j) T = C^cf), where 

c ~=-h4 +l - (12) 

and so, in the space of polynomials, we have now 2k expanding directions corresponding 
to £ n , for n < 2k. 

Equation (11) is solved by putting 0(£, r) = e~p irj h(e~ T ^ 2k ^,T) whereby 4-h{y : T) = 
e~ T h(y, r) p and so, for p = r — r 

(P{ ^ ,T> [l-(p- l)f(e-P/ 2k ^)P- x {l - e-p)] 1 ^- 1 { ' 

where (f>(£, r ) = /(^). The stationary solutions / fe * are stable in a suitable codimension 
2k space: let us consider / smooth, with 

/(0) = {p- l)-^i, /W(0) =0 £ < 2fc, /^(0) = (3 < (14) 

and 

o</(o<(p-i)"^ e^o. (is) 



Then, for all (eR, 



where 



m,r)-m)\-io (i6) 



T^OO 



A*(0 = (p-i + &r fe )"— (17) 

for some b depending on fi,k, p. 

These considerations thus lead us to expect (8) to have global solutions with initial 
data in a suitable codimension 2k set in a ball around (17) in a suitable Banach space. 
Of course, the perturbation L~ 2 (p^ in (8) is very singular, but, basically, the picture is 
not much modified: the unstable modes turn out to be r-dependent Hermite functions 
instead of the monomials £ n , and one has to fine-tune the projection of the initial data 
on these Hermite functions instead of imposing the vanishing of the derivatives as in 



(14). See |15J for details. 



4.3 A non- conventional center manifold problem: k — 1. 

Consider now the case k = 1. There are several differences with respect to the previ- 
ous one: there is no damping factor in front of the second derivative in (8), and the 
asymptotics is given by (4), with a "universal" b*, and an extra logarithmic factor in 
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the definition (3) of g(t,T). To understand the dynamics of (8) in that case, let us start 
by considering again its linearization around the constant solution <fi = (p — 1)~. The 
linear problem is <p T = £<p, where now 

Hence, the first thing we have to do, in order to understand the stability of the constant 
solution, is to study the spectrum of the linear operator C. C is self-adjoint on T>{£) C 
L 2 (R,dfi) with 

dKO = —^ (19) 

V47T 

The spectrum of C is 

Ti 

spec(C) = {1 - - | n E N} (20) 

and we take as eigenfunctions multiples of Hermite polynomials 

mo = E • (-lrr 2 " (21) 



that satisfy 



and 



"n"m"'IJ' ^ Tl-Onm \^"^) 



Ti 

Ch n = (1 - -)h n . (23) 

Thus the linearization of (8) at the constant solution, for k — 1, has two expanding 
("relevant") modes, ho and hi, and one neutral ("marginal") one, h 2 = C, 2 — 2. 

Our goal is therefore to construct a center manifold for the flow of (8), in a neigh- 
bourhood of the fixed point. Formally, we would expand, 

oo 

0(f , r) = (p - l)^r + ^ T ) as $ = J2 Mr)hn(Z), (24) 

n=0 

and rewrite (8) for the ip n {T~) as an infinite set of ODE's: 

—i^nij) = (1 — —)ip n (j) + nonlinear terms 
dr 2 



A formal solution of this flow yields (see below for the calculation): ip2(j) — C p (\og t)t l 
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with C p = -b-V p as in (4). 



However, there are severe problems with this approach. Since the eigenfunction h n of 
the linearization C increase at infinity, the expansion (24) is not useful for £ large, and, 
in particular, we cannot use any standard infinite dimensional center manifold theorem. 
The key to the solution to this problem comes again from a scaling argument, which will 
explain the emergence of the fixed point /*. Let 

L (£,r) = 0(L£,L 2 r) (25) 

Then, (pi satisfies the equation 

if(0 L ) = L- 2 (-0 Lr + L ^), 

where we defined 

H{<f>) = \^ + ^-<p-<P P ■ (26) 

2 p — 1 

Hence, as L — ► oo, we expect the solutions of H{cf) = to be relevant. These are (like 
the stationary solutions of (11)) given by the one-parameter family f£, given by (6), 
with k — 1. Therefore, we have the following picture: instead of perturbing around the 
constant solution, introduce </>&(£, r) = (p — 1 + &£ 2 / r )~ , and write 

<f>(Z,T)=MZ,r)+ri(Z,T). (27) 

A local (i.e. for £ small) center manifold analysis then fixes b = b*, as in (4). For large £, 
namely — > 0(1), the linearization C of (8) around 0& differs from C, and, actually, in 
that region (see below) C ~ £ — 2. By (20), the spectrum of C — 2 is entirely negative. 
Hence, the dynamics tends to contract rj. 

We will now explain the calculation that yields b* in (4). For simplicity of notations, 
we shall consider only p = 2, i.e. 



We get, using H(tp b ) = 0, 



fe (£,r) = (l + ^ 2 /r)- 1 . (28) 



r}^ - H((j) b + i]) + H(<f> b ) + b ^ - (j) bT 

(C + W) V + M( V ) + (j)^ - (j> bT (29) 



where we introduce 



W = 2(0 6 - 1) (30) 

M( V ) = ( ( j )b + r ] ) 2 - ( j ) 2 b -2<t )b r ] . (31) 

The operator £., given by (18), has two unstable modes. Note that, formally, (i.e., for £ 
of order one) W is (9(r _1 ), M is nonlinear in rj and (f) b ^ — <p bT is 0(t~ v ). We want to 
construct a center manifold for (29), i.e. to see how to fix two parameters in the initial 
data of rj, such that the flow of (29) drives i] to zero. A simple calculation will show 
that this can only be achieved through a suitable choice of b in (28). 
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To see this, let us simplify further and consider 77 even in £, which will imply that 
we need to fix only one parameter. This example contains all the relevant features of 
the general case. It is convenient to write 

Vo(r) = - (32) 

T 

and define i[) by 

77 = i] + if). 

Then ip satisfies the equation 

i> T = Ci) + N(i)) + a (33) 

with 

£ = c + v 

V = 2(0 b + ?7o-l) (34) 

N(iP) = (0 6 + r ?o + V) 2 -(0 fe + %) 2 -2(06 + %)^ (35) 

a = (fr b {:{: - (p bT + (£ + W)i] - r] 0r + M(r] ) 

= Mi -'Pbr + Vo + Wrjo -r] 0r + M(r] ). (36) 

Let us decompose ip (which is even in 77) as 

il> = MT)+MT)h 2 + il> ± (37) 

where ip 1 - is orthogonal to h n , n < 2 . Next we expand V and a (for £ = 0(1)): 

2be 2 2a „, 1 N 
V = ^- + — + 0(— ) (38) 

T T T z 

a = (a - 26)r" 1 + (a + a 2 + (12& 2 - 6 - 2ab)£ 2 ))r- 2 + C(r- 3 ). (39) 

Inserting (32), (37) in (33) and retaining only the leading terms in 1/r and ipi, i = 0, 2, 
we get from ipi T = (2 l iX)~ 1 {hi,ip T ) ((•,•) is the scalar product of L 2 (R,dfi) and we use 

(22)): 

</>0r = ^0 + (a - 2b)r- 1 + i? (40) 

^ 2 r = /3t-V 2 + (126 2 - b - 2ab)r' 2 + R 2 (41) 

where i? = 0{t~ 2 + r" 1 ^! + |V| 2 ),i? 2 = C(r~ 3 + r" 1 !^! + t- 2 |^ 2 | + |^| 2 ),and /3 = 
2a — !&(£ 2 /i2, /12) = 2a — 20b (coming from the Vip term in (33)). We choose now a so 
that the 0(r~ l ) term in ipo T vanishes i.e. 

a = 2b (42) 

and b such that the 0(t~ 2 ) term in ip2 T is zero: 

b = b* = 1/8. (43) 
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Note that this choice correspond to b = b* in (4) for p = 2. Then (3 = —2 and our 
equations read 

2 

IpOr = *Po + Ro, V>2r = 02 + #2- (44) 

T 

Now, keeping in mind the presence of the Rq, R<i terms, 

*P = (9(r- 2 ), V>2 = G((log r)r- 2 ) (45) 

would be consistent solutions. Of course, we need to show that the expanding variable 
■00 wm satisfy (45) by a suitable choice of ipo{r ). This is rather easy to do, using the 
fact that ipo is expanding; in the general case (with more than one parameter to fix), we 
used a topological argument. 

In the rigorous proof, we set up a suitable Banach space for the function rj, parametrized 
by the ipi's for £ 2 /r < 0(1), and a function r\i for £ 2 /r > 0(1). The function rji will 

contract under the action of C: indeed, from (34,28,32), we see that the potential V 
tends to —2 as £ 2 /t (and r) — > oo. 

5 Open Problems 

There are several patterns and fronts in dissipative equations, and their stability can 
probably be studied using RG methods. For example, the Cahn-Hillard equation [|21"|j : 



u t = A(-Am - u - u 3 ) (1) 

is often used to study the phase separation in alloys and fluids. One would like to study 
the stability (and possibly the dynamics), in infinite volume, of interface solutions of 
that equation. 

Another major open problem consists in the extension of the RG method to hyper- 
bolic equations. The latter have their own scaling laws, see [^, |6^] . A lot is known on 



the stability of soliton solutions of (generalized) KdV equations [5B|, but the stability of 
localized solutions of other nonparabolic equations is quite open. Also, there are open 
problem concerning the blow-up of solutions, most notably in the nonlinear Schrodinger 
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So far, we have only considered equations whose solutions have a rather simple 
asymptotic behaviour. Of course, it is well known that finite dimensional dynamical 
systems described by differential equations can have a chaotic asymptotic behaviour, 
i.e. depend sensitively on the initial data but have also some statistical regularity in 
the sense that the long time average along the orbits is described by an invariant (SRB) 



measure |27|]. For certain classes of F in (1.1), one would like to find natural invariant 
measures for the flow. 

A class of dynamical systems, modeling PDE's, is obtained by discretizing space and 
time and considering a recursion 

u(x,t + l) = F(x,u(-,t)) (2) 
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i.e. u(x,t + 1), with x being a site of a lattice, is determined by the values taken by 
u at time t (usually on the sites in a neighbourhood of x). For a suitable class of F's 
such dynamical systems are called Coupled Map Lattices (CML) [EI], [52J. The study of 



the invariant measures, and their properties, even for these "toy models" is essentially a 
terra incognita. In |2U], |5D, [T7| various "high-temperatures" (weak coupling) results are 



obtained. Almost nothing is known about the "low temperatures" (strong coupling) side 
iP|, 1571). One would like to find models for which the set of invariant measures 



[see | 

changes as the coupling is varied; this phenomenon would then be interpreted as a kind 

of nonequilibrium phase transition. 
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